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We present a theoretical and experimental investigation of the emission characteristics and the 
flux of photon pairs generated by spontaneous parametric downconversion in quasi-phase matched 
bulk crystals for the use in quantum communication sources. We show that, by careful design, 
one can attain well defined modes close to the fundamental mode of optical fibers and obtain high 
coupling efficiencies also for bulk crystals, these being more easily aligned than crystal waveguides. 
We distinguish between singles coupling, 7s and 71, conditional coincidence, and pair coupling, 
7c, and show how each of these parameters can be maximized by varying the focusing of the pump 
mode and the fiber-matched modes using standard optical elements. Specifically we analyze a 
periodically poled KTP-crystal pumped by a 532 nm laser creating photon pairs at 810 nm and 
1550 nm. Numerical calculations lead to coupling efficiencies above 93% at optimal focusing, which 
is found by the geometrical relation L/zr to be « 1 to 2 for the pump mode and « 2 to 3 for the 
fiber-modes, where L is the crystal length and zu is the Rayleigh-range of the mode-profile. These 
results are independent on L. By showing that the single-mode bandwidth decreases oc 1/L, we 
can therefore design the source to produce and couple narrow bandwidth photon pairs well into 
the fibers. Smaller bandwidth means both less chromatic dispersion for long propagation distances 
in fibers, and that telecom Bragg gratings can be utilized to compensate for broadened photon 
packets — a vital problem for time-multiplexed qubits. Longer crystals also yield an increase in fiber 
photon flux oc and so, assuming correct focusing, we can only see advantages using long crystals. 

PACS numbers: 03.67.Mn, 03.67.Hk, 42.50.Dv, 42.65.Lm 



I. INTRODUCTION 

Spontaneous parametric downconversion (SPDC) ac- 
counts for the majority of entangled photon pairs be- 
ing produced today. It can be described as a process 
in which the electromagnetic field of a single photon — 
traveling inside a dielectric material such as a bircfringcnt 
crystal — interacts with the atoms by absorption and gives 
rise to a nonlinear response in the field of polarization, 
thereby leaving the possibility of two or more photons be- 
ing re-emitted. The laws of conservation of energy and 
momentum, together with the randomness and indistin- 
guishability in the process, also give rise to entanglement, 
a nonlocal correlation between the photons. 

In quantum communication numerous experiments 
have been performed to date involving non-cntanglcd or 
entangled photons being sent over longdistances, e.g., 
sources of heralded single photons P, HjBi quantum 
cryptography ^ |^ , and teleportation [J- A typical 
such experiment involves launching each photon of a (en- 
tangled) pair into single-mode fibers and to deliver each 
one to a separate party for encoding or decoding. For 
successful distribution over long distances it is vital to 
have a high rate of pairs generated at the source, as the 
attenuation of the fiber is a strongly limiting factor even 
at the wavelength of 1550 nm for which the fiber is most 
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transparent. Today, results with crystals of periodically 
poled materials have proved this viable even at moder- 
ate pump laser powers Q , and in some cases the problem 
has turned into a matter of limiting the pump power to 
avoid creating two pairs at the same time, as this will give 
false coincidences also when having low smgZe-coupling 
efficiencies. Instead, what has gained importance is to 
have a high pair-coupling efficiency that increases the 
probability of both photons of a pair being present in 
the fibers once they have been created. Furthermore, the 
use of time-multiplexed schemes |^ ^| have elicited the 
need of launching photons having very narrow frequency 
bandwidth and long coherence length in order to limit the 
effects of dispersion in the fibers, and to enable the use 
of interferometers. Rather than just filtering the emis- 
sion at some desired width, as is commonly done, we will 
show that it is more efficient in terms of photon-rates to 
design the source so that the bandwidth is determined 
by the crystal length and fiber coupling alone. 

It is the purpose of this article to calculate the maxi- 
mum coupling efficiency achievable for photon pairs gen- 
erated in crystals that are phase-matched for colinear 
emission in general, and for periodically poled KTiOP04 
(PPKTP) crystals using non-degenerate quasi-phase 
matching (QPM) in particular. Wc look for the opti- 
mal condition for focusing of the pump onto the crystal 
and focusing of the emission onto the fiber-end (mode- 
matching) which maximizes cither the single or the pair- 
coupling efficiency. The focusing is specified using the pa- 
rameter ^ — L/zr, adopted from [Il| with a slight mod- 
ification, where L is the length of the crystal and zj^ is 
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the Raylcigh range. We make no thin-crystal approxima- 
tions, but take fully into account the focusing geometry 
of all three interacting fields: pump, signal, and idler, by 
decomposing all three fields into a complete set of orthog- 
onal plane- wave modes. Other optimizable parameters of 
these beams include the direction of the beam axis and 
the location of the focus. Both are regarded fixed, the 
former being motivated by the colinear geometry of per- 
fect quasi-phase matching, and the latter by the fact that 
focusing onto the center of the crystal shows to give high- 
est efficiencies. (Support for the last claim is given in [T3I 
for second harmonic generation.) We also regard the cen- 
ter frequency of the beams, the power of the pump, and 
the optical properties of the crystal as fixed parameters 
of the problem. We take into account the polychromatic 
character of the emission but assume a monochromatic 
pump (continuous- wave pump), and we investigate how 
the coupling efficiency depends on the length of the crys- 
tal and the bandwidth of the wavelength filter in front of 
the fiber, but also how the fiber coupling affects the band- 
width of the coupled photons and the achievable photon- 
rates. Our goal is to give a simple recipe for setting up a 
colinear source of entangled photon pairs that optimizes 
the focusing for the highest single and pair coupling effi- 
ciencies into single mode fibers, and that also determines 
a suitable crystal length for a desired bandwidth. 

Shortly after the demonstration of parametric gener- 
ation (PG) and second harmonic generation (SHG) in 
the 1960s, Boyd and Kleinman [ij, and others, ad- 
dressed the focusing in non-colinear geometries of type-I 
and showed the importance of optimization for achiev- 
ing maximal conversion efficiency in optical parametric 
oscillators and frequency doublers. By using cavities to 
enhance the processes one can control the spatial mode 
of the pump, signal and idler to support only the fun- 
damental TEMqo mode, and under this condition Boyd 
and Kleinman suggested that the general optimal fo- 
cusing is to set the ^-parameters of all fields the same 
(^p = = 6)- Later, Guha et al. ^3 showed that hav- 
ing unequal parameters can improve the conversion even 
further and this is also supported by our results. The 
case of type-II SHG have also been studied 0|, as well 
as sum- and difference frequency generation (SFG and 
DFG) 111, with similar results. These works were all 
treating the light as a classical field, having the signal 
beam acting as the relatively strong control-field that 
is being amplified by the much stronger pump-beam to- 
gether with the creation of an idler. It is not unreasonable 
to expect that a different situation arises at the quantum 
level where both the signal and idler initially are in un- 
controlled vacuum-states. 

Spontaneous parametric downconversion commonly 
takes place in bulk crystal configurations where the sig- 
nal and idler modes are not restricted by cavities. This 
will provide an additional degree of freedom. The pump 
is assumed to be TEMpo, but the emission will in gen- 
eral be spatially multimode. A central problem in this 
article is to find how much of the emission is in a trans- 



verse and longitudinal fundamental singlc-modc at differ- 
ent focusing conditions. For the transverse part, such a 
single-mode, being Gaussian shaped, is very close to the 
Bessel function of the first kind, Jo (a), which describes 
the shape of the fundamental fiber mode, and will there- 
fore provide nearly perfect overlap. After determining 
the mode of the emission we also calculate the fac- 
tor, commonly used as a measure of beam-quality, and 
compare it to experimentally obtained results. 

To our knowledge, no analysis has been made to date 
that characterizes the colinear emission in quasi-phase 
matched materials in the way presented here, i.e., mak- 
ing no assumptions about short crystals or weak focusing. 
It should be noted that the analytical calculations be- 
come difficult without these assumptions and so our goal 
have been to formulate the final expression in such a way 
that it can be evaluated numerically with relative ease, 
with only simple assumptions being made. Taking into 
account all the needed degrees of freedom — azimuthal 
and polar angular spectrum and frequency included — 
these numerical computations will become quite time- 
consuming on an ordinary personal computer, but still 
doable. 

Various other attempts have been made in the past 
to characterize the one- and two-photon spatial opti- 
cal modes generated by non-colinear birefringent phase- 
matching. However, most of them do not use single-mode 
fibers to collect photons; Monken et al. and Pittman 
et al. ^3 show how focusing of the pump with a lens can 
increase the coincidence counts using an analysis limited 
to thin crystals, and Aichele et al. [13 seek to match the 
spatio-temporal mode of a conditionally prepared photon 
to a classical wave by spectrally and spatially filtering the 
trigger, however, without considering focusing effects. 

More recent work connected to ours is a number of pa- 
pers that consider the coupling into single-mode fibers; 
Kurtsiefer et al. provide, for thin crystals, a hands- 
on method of determining the mode of the emission 
using the relation between the emission-angle and the 
wavelength coming from the phase-matching conditions. 
For maximal overlap between the emission-mode and the 
fiber-matched mode (target) they presume it is best to 
choose the waist of the pump-mode and fiber-matched 
mode equal. According to [ll|, and our results, this is 
not optimal in general. Bovino et al. |20l | take on a more 
sophisticated approach as they carry out the biphoton- 
state calculation for a non-colinear source, which takes 
into account focusing, dispersion, and walk-off and ar- 
rives at a closed expression for the coincidence efficiency. 
Other work have been continued along the same fines 12 ll : 
our conclusion from examining the formulas herein being 
that high efficiency can always be achieved for any length 
of crystal by choosing the pump waist large enough and 
the fiber-matched waist small enough. This is in contrary 
to our results which show an optimal value of the focus- 
ing parameter (1 ^ ^ < 3 ). Furthermore, as shown both 
in this report and in |lll |. for a specific crystal type and 
wavelength configuration the value of ^ is found to be 
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a fixed constant for all crystal lengths which makes the 
pump-beam waist Wq relate to the length as W n oc VL 
(at optimal focusing), while the results of Ref. ]2nL I2H 
appear to show a linear relationship. We are not sure 
whether these apparent differences are best explained 
by the different situations of a non-colincar and colinear 
source, pulsed vs. continuous-wave pump, or by other- 
wise different models or parameters in cither case. It can 
be noted that our results seem to provide good agreement 
with experiments. 

The particular source of photon pairs that spurred the 
work of this article is presented by Pelton et al. in Ref. 
[2^ . The main idea is to create polarization-entangled 
photon pairs at the non-degenerate wavelengths of 810 
nm and 1550 nm from a purnp-photon at 532 nm, us- 
ing two orthogonally oriented [23, long, bulk KTP crys- 
tals. These crystals are periodically poled for quasi-phase 
matching which provides colinear emission suitable for 
coupling into single-mode fibers, but as told, also require 
some optimization for maximum throughput. Prelimi- 
nary results can be found in . Related work is found 

in lElMEi 

The agenda of this article is as follows. Section |n] 
gives a mathematical background, starting in subsection 
III Al with a review of the one and two-photon state of 
the emission derived in Appendix IXI In subsection III Bl 
we calculate the emitted modes, which arc qualitatively 
measured using the beam quality parameter JVP . This is 
followed in subsection III CI bv a mathematical definition 
of the single-coupling, coincidence, and pair-coupling ef- 
ficiencies. Section IIIII presents the numerical results of 
the couplin g (|III Al - ITTTbI) . bandwidth IjlllCp . and the 
AP factor Ijlll pp . Section Hvl covers the experimental 
setup and the experimental results, where a comparison 
is made to numerical predictions. Conclusions are found 
m Section El 



II. THEORETICAL DESCRIPTION 

The aim of this section is to derive the formulas used 
for the numerical calculations of the emission modes, cou- 
pling efficiencies, and emission bandwidths for the emit- 
ted quantum state of the SPDC process, and also to give 
a physical meaning to these concepts in the role of single 
photon sources. We will optimize over the spatial param- 
eters involved to find the highest quality modes and max- 
imal coupling efficiencies attainable. The result is based 
on a calculation carried out in Appendix fXl involving the 
Hamiltonian that governs the interaction of spontaneous 
parametric downconvcrsion in quasi-phase matched ma- 
terials. The crystal is pumped by monochromatic and 
continuous wave laser light (p) of frequency Wp, which is 
propagating in a Gaussian TEMqo mode along the z-axis, 
producing a signal (s) and idler (i) field in the same direc- 
tion. Figurendefines the laboratory axes used; the z-axis 
being along the length L of the crystal, the x-axis along 
the height, and the y-axis along the width. The crystal 
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FIG. 1: The figure sliows the periodically poled crystal with 
the laboratory coordinate system drawn. Also defined are the 
crystal's axes X, Y, and Z, referring to the polarization of the 
incoming and outgoing electromagnetic fields. 



is bi-axial, and the crystal axes X, Y, and Z are ori- 
ented as shown in the figure. We have chosen the poling 
period in the crystal to allow for co-polarized (ZpZgZi), 
colinear down-conversion, but the calculations are gen- 
eral enough to allow other polarization settings. The 
refractive indices, and thus the phase-matching, depends 
on the temperature of the crystal and is determined by 
the Sellmeier coefficients of PPKTP [Ullll. In general 
we are interested in phase-matching at non-degenerate 
wavelengths, and for such cases the shorter wavelength 
will be regarded as the signal and the longer wavelength 
as the idler. 

Many references, following Klyshko 12^ , start with the 
coupled mode equations and look at the evolution of op- 
erators to find the two-photon state from SPDC in terms 
of a frequency and angular intensity distribution [30| . 
This is effectively the same as finding the diagonal ele- 
ments of the second order moment density matrix which 
represent the incoherent part of the information of the 
state. This information is sufficient for determining the 
shape of the emission. However, it is not sufficient for de- 
termining the overlap between the emission and a single- 
mode fiber. In this case we need the "coherent" informa- 
tion available in the full density matrix. The approach 
we take in Appendix 1X1 and in the next subsection is to 
use the Schrodinger picture and look at evolution of the 
state to find the two-photon amplitude. In the following 
subsections we then diagonalize the corresponding den- 
sity matrix into a sum of coherent parts (eigenmodes), 
and project each one onto the fiber-mode so that wc can 
calculate the coupling efficiency as a sum of overlap coef- 
ficients. We also use this decomposition to calculate the 
electrical field and beam profile of the emission. 

A. The emitted two-photon state 

The two-photon amplitude describes the joint state of 
the signal and idler emission in terms of (internal) an- 
gular and frequency spectrum. Using spherical coordi- 
nates (see Fig. ^| the two-photon amplitude derived in 
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Eq. HA28|I becomes 



47r2x2/iL 



ih 
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X sine 



(1) 



where, according to Eq. I|A26(I 

Afc; = fcs cos 9s + h cos 6*1 - ^/l - {P^ + Q^) + /\ , 

(2) 

and, according to Eq. HA25|) 

p2 + = 

fcg sin^ 0s + fcf sin^ 0i + 2fcsfci sin 9s sin 0i cos(A(p) 



{k^Y 



(3) 



All three interacting fields have been decomposed into 
a complete set of orthogonal plane-wave modes, k{9,(p). 
The magnitudes of the /c- vectors, kg and fcj, arc given 
by Eq. (|A17p . 9s and 6*1 are the internal polar angles 
of the plane waves of signal and idler respectively, Aip 
is the difference in angle between the azimuthal angles 
ifs and ifi, and e is the frequency (specified by a single 
parameter due to exact energy-matching). Furthermore, 
X2 is the nonlinear coefficient of the crystal, K is the 
grating constant of the poling, L is the length of the 
crystal, and Wq^ is the pump-beam waist radius. A(e) is 
the frequency amplitude of the detector filter having a 
bandwidth AA (FWHM) and a center wavelength (all 
wavelengths in vacuum). Via the relation e = 2TTc(n^/ X — 
/Aj,) its form, assuming a Gaussian shaped filter, is 
given by 



Aie;X) 



-21og(2)(A-A^)^/AA^ 



(4) 



In a plane wave mode-decomposition, Eq. ^ repre- 
sents the two-photon field (that is generated in the crys- 
tal by the pump field) in the form of a continuous an- 
gular spectrum in polar and azimuthal degrees of free- 
dom. Together with the frequency, the full state is a 
tensor-product of four degrees of freedom. We will need 
to discretize the spectrum in order to represent it on a 
computer. As the size of the Hilbert space of the full ket- 
vector becomes very large for a large number of points in 
resolution, we need to limit its size to make the numeri- 
cal calculations feasible. In the following, the two-photon 
state is therefore explicitly represented only by the polar 
angles of the signal, 16*3), and the idler, 16*1), written as 
kets, leaving the state implicitly dependent upon the two 
remaining degrees of freedom, Aip and e. The purpose of 
this notation is to refiect the actual way that the state is 
numerically implemented as a one-dimcnsional array of 
9 (the density matrix is a two-dimensional array), with 



separate arrays being calculated for each discrete value 
Aip and e. Choosing Ng discrete plane-wave modes as a 
basis of the polar angle, the two-photon state can then 
be formulated as 
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(5) 



There are a few approximations that have been made 
during the calculation of S, apart from the paraxial ap- 
proximation inherent in the standard form of the angular 
spectrum representation of the Gaussian pump field of 
Eq. ljA19|l . These include: i) the assumption of a con- 
stant pump k-vector magnitude kp = k^ in order to re- 
move the implicit dependence of 9p and (pp in Eq. ljA16p . 
which thus leads to Eq. ii) the assumption of an 
infinite coherence length of the pump (cw), providing a 
(5-function over frequency so that we can describe the 
signal and idler by a single frequency e, and Hi) the as- 
sumption of having the same refractive indices along the 
crystal's X and Y axis, such that the X-component of 
the A:- vectors can be set to the same as that of Y. The 
last assumption also provides a motivation for the out- 
put of completely rotationally symmetric modes, and will 
greatly simplify the expressions and the numerical calcu- 
lations as the azimuthal angle dependence, via ips and 
(pi, is automatically removed from the two-photon ampli- 
tude. The two-photon density matrix is given by 



(6) 



which now contains four degrees of freedom; 9s and 9i 
being the two state parameters, and A(p, e being two 
other parameters which we will trace over later. Note 
that psi is a description of the emission inside the crystal, 
not taking into account the refraction between crystal 
and air. 



B. The emission modes and the beam quality, M 

We are interested in the shape of the signal or idler 
beam profiles using free detection so that we can com- 
pare with images taken by a CCD camera. To do this 
comparison we need to have the beam described in terms 
of the electrical field, which is given as the Fourier trans- 
form of the angular spectrum (the density matrix). The 
electrical field, or intensity, then gives the beam profile 
which, in turn, determines the factor. 

First, each signal or idler beam are made independent 
of the other beam by partially tracing over its partner. In 
the following we trace over the signal in the polar angle 
degree of freedom, and in doing so we get the reduced 
density matrix for the idler. 



Pi 



Aip.e 



Ng 



(7) 
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The remaining dependence on A<p can also be removed 
following the standard trace-operation, which is here 
equivalent to a sum over density matrices, 



Emission modes 



Pi 



TrAv(pf'^'') = 



(8) 



Additionally, as we could in principle measure the fre- 
quency of the photons at a resolution given by A Arcs = 
A^/cAtgate (set by the timing information of the detec- 
tors, > 1 ns, to be < 8 pm), which generally is much 
smaller than the bandwidths of the filters, we need to 
incoherently sum over the frequency e in the same way, 
giving a final pi describing the state of the idler. 



Pi = Tre(pf) = 



(9) 



Mode decomposition 

We cannot, however, directly now apply a Fourier 
transform to the reduced density matrix pi, as it is gen- 
erally mixed. Instead, we shall diagonalize pi to find 
its eigenvectors and eigenvalues. For such a Hermitian 
matrix all eigenvalues are real and the eigenvectors will 
form a complete orthonormal set. Thence, the set will 
represent a natural mode-decomposition of the emission, 
and consequently, each vector, or mode, will represent 
a coherent part of the emission. The sum of all modes 
weighted by its corresponding eigenvalue will determine 
the state. For each such mode, on the other hand, we 
can apply a Fourier transform and thus find the electri- 
cal field modes. The squared sum of all electrical field 
modes, again weighed by the corresponding eigenvalue, 
will then determine the total electrical field. We will 
quantify this to show our future notation; the reduced 
density matrix is first diagonalized by T'^pT = D, such 
that T ~ (ICi)) ^2)7 • ■ • J ICa^o)) has the eigenvectors in 
the columns, and D has the eigenvalues A„ in its diago- 
nal elements. The result is a density matrix that can be 
represented as a sum of pure states. 



Ng 

P = ^ A„|C«)(Cn 



(10) 



where Ng is the Hilbert-space dimension. Following this 
result, in Fig. [3 is plotted the one-dimensional angu- 
lar spectral form u[0y], taken as an integration of the 
absolute square of the two-dimensional angular spectral 
amplitude a2,y[6']. We have a^ylO] — J2n ^riCn[()], where 
Cn[0] is the discrete function representation of \Cn), and 
0^=el + el. Hence, 
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y 



(11) 




FIG. 2: The figure shows an example of the angular spectral 
form u[ey] of the emitted idler light att 1550 nm in a PPKTP 
crystal (central curve) which gives an factor less than 3 
with a filter bandwidth AA = 10 nm. The pump at 532 nm is 
focused close to optimal, = 1.3. The insets show the four 
lowest order LGpo modes which are similar, but never the 
same as the natural eigenmodes of the emission, and illus- 
trates how the factor in general grows with mode order. 



The field intensity 

We can now transform the angular spectrum modes 
\(n), into electrical field modes En- As these modes 
are rotationally symmetric and depend on one param- 
eter only, the electrical field is most suitably expressed 
through the Hankel transform. In writing the trans- 
form in the following form we make use of the fact that 
the vector iCn), again written as a discrete function, 
Cn[^7 v] = Crjfi*]; is independent of ip. Thus, 



~ikz cos 6 



Jo (^ky/a 



(12) 



where the basis functions Jo{a) of the Hankel transform 
are the Bessel function of zero order and the solution to 
^ J^^ exp {i a cos (fi) dip. However, the one-dimensional 
Fast Hankel Transform (FHT), which would possibly pro- 
vide very fast computations, is not widely implemented, 
at least not in an efficient form for use in Matlab or 
Mathematica and was not available to us at the time 
for the numerical calculations. Therefore, the next sim- 
plest transform at hand is the two-dimensional Fourier 
transform. 



e V 



-ikz cos 6 



is the one-dimcnsional angular spectral form. 



^ ^kxsinOcos^p ^kysinOsin^p (13) 



6 



With still two dimensions being used, Eq. (|13|l can also 
be rewritten using the polar angle components 9^ and Oy^ 



En{x,y,z) = ^^KCn 



X e e " 



g-»fcz cos(^i|Tef ) 
(14) 



where 



+ 6*2 . In this form, which is the form 

we will use, Eq. p4|l represents a standard single two- 
dimensional FFT. Note that this transform is, in general, 
not separable with respect to x and y into two, but sim- 
ple, one-dimensional transforms. This is a characteristic 
of Laguerre-Gaussian modes and of the modes emitted 
by the crystal, in comparison to Hermite-Gaussian modes 
which are always separable. 

The intensity is now given by incoherently summing 
all field-modes. 



I{x,y,z) = ^ \En{x,y, : 



(15) 



Finally, the transversely integrated intensity profile of the 
emitted beam is given by /(y, z) = I(x, y, z). 



Gaussian beam fitting 

The beam waist radius w{z) can be found from the 
standard deviation a{z), or the second moment, of the 
intensity distribution I{y,z), as w{z) ~ 2(t{z), see Ref. 
|3l| . The standard deviation is known to provide the 
correct waist estimate for arbitrary multimode light as 
opposed to trying to make a curve-fit with various mode- 
shapes. Readily, cr^(z) = T,y {v - y{z)f Hv, z), where 
y{z) = J2y uHy^ ^) is the expectation value with respect 
to the spatial position y in the intensity distribution. As 
said, we will use the beam quality factor to quan- 
tify the emission. This factor is determined through the 
Rayleigh range 



entering the standard Gaussian beam formula 



ZR 



(16) 



(17) 



By varying the parameters and we can make a 
curve-fitting of the model profile Wmodoi(-z) to the actual 
beam profile w(z), such that the factor is determined. 
Equation states that the diffraction limited funda- 
mental Gaussian mode TEMqo has a beam quality factor 
of = 1. As a comparison, this factor increases for gen- 
eral higher order Laguerre-Gaussian modes LGpm |32l |. 
defined by the radial index p and the azimuthal mode 
index m = 0, such that = 3 for p = 1, AI^ = 5 for 
p = 2, and = 7 for p = 3 and so on, see Fig. [3 
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FIG. 3: The figure shows a Venn diagram. It illustrates the 
single coupling efficiencies 7s and 7i, pair coupling 7c, and 
conditional coincidences and which are defined in 
the text. The total amount of pairs Q,p generated within the 
bandwidth of the detector filter AA is normalized to unity, 
and represents perfect coupling. 



C. Single coupling, coincidence, and pair coupling 

To characterize the source and to optimize the cou- 
pling of the emission into optical fibers we shall make 
use of three parameters: single coupling, conditional co- 
incidence, and pair coupling. However, before we define 
each of the three coupling parameters we shall briefly 
comment on the necessity to relate them to the detection 
window being used, i.e., the frequency bandwidth of the 
detector filter AA. The emission will always fluoresce in 
a wide spectrum, and in that sense there is no mean- 
ing to speak about a coupling efficiency for photons that 
cannot be seen through the window in any case. By mak- 
ing a simple normalization to the filter bandwidth, the 
coupling probability will consistently measure only how 
well photons of specific frequencies are spatially collected 
into the fibers. For example, for any fixed filter and no 
spatial filtering, as is almost the case with a multimode 
fiber, and certainly the case in free-space, the coupling is 
always perfect. Effectively, this normalization enters the 
calculations through the bandwidth in Eq. Q). Figure |3| 
helps to illustrate the different coupling parameters using 
a Venn diagram. 



Single coupling 

The single- coupling efficiencies 7s and 7i are readily 
defined as the probability to find a photon in the fiber 
which has been emitted within a certain filter bandwidth. 
The single-coupling efficiency is useful when maximizing 
the individual rate of photons present in the fibers. To 
calculate the probability we shall take the overlap of the 
emitted modes with the mode of the fiber as seen from 
the crystal, here called the fiber-matched mode. That is 
to say, the form of the mode that can be traced back to 
the crystal from the fiber-tip, not worrying about crys- 
tal refraction or any other optics in between performing 
the actual transformation. Also, we do not consider any 
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FIG. 4: The picture shows the geometry of focusing, with 
the Rayleigh-range zr, the crystal length L, the beam waist 
radius wo, and the focus offset zo being defined. The focusing 
parameter is defined as ^ = L/zr. 



^°P*). In both Eq. O and Eq. ^ the parameter ^ 
enters through the beam waist radius of the pump mode 
wop and the signal/idler fiber-matched mode wqq, accord- 
ing to wop = ^/iAp/TT^p, and wqo = ^jLX~pK^^. We 
can formalize the optimization of the signal and idler 
fiber-matched modes as 

7°P' = inax 7(ep,Cs.i), (20a) 

SS,i 

rP* = argmax 7(ep,4.i), (20b) 
with 7 given by Eq. H19|). 



additional aperture limitations enforced, e.g., by irises. 

The true mode of the fiber is described by a Bessel 
function. However, it can be approximated very well 
with a fundamental Gaussian which in normalized form 
is described by 

•\/27r 

where Wqq is the beam waist radius of the fiber-matched 
mode, TEMqo, as determined by the focusing system, 
and Zqq is the location of the corresponding focus (which 
shall be at the center of the crystal Zqq = for optimum 
couphng), see Fig. 0] 

The single-coupling efficiency is trivially given by 7 = 
Tr(|G'oo)(Goo|p), but the numerical optimization con- 
verges slowly and badly using this form. For this reason 
we shall exploit the diagonalization and calculate the sin- 
gle coupling efficiency as the sum of the projection of each 
emitted mode |C„) onto the fiber-matched mode |Goo), 

71 = 1 

where jCn) is given by the density matrix, ps or pi, as 
defined by Eq. (|10|l . resulting in 7s or 71 respectively. 

Optimization 

The maximum achievable coupling efficiency is deter- 
mined by an optimization of Eq. H19(l with respect to 
the focusing conditions of either the pump mode, or the 
fiber-matched signal/idler mode, or both. To quantify 
the focusing we shall use the beam focusing parameter 
^ = L/zn, where L is the length of the crystal and 
zb is the Rayleigh-range (note that we have = 1 for 
both the pump mode and the fiber-matched modes). See 
Fig. ^ The parameter is suitable as a dimensionless rep- 
resentation of the focusing geometry. (As will be shown 
further ahead, the results indeed show that the geome- 
try is kept intact at optimal focusing, irrespectively of 
the length of the crystal, which corresponds to a fixed 



Conditional coincidence 

The conditional coincidences, fis\i and ^i|s are useful for 
the characterization of heralded single photon sources, 
and are defined as the probability to find a photon in 
either the signal or the idler fiber given that the partner 
photon has entered its fiber, whether or not its detected. 
The conditional coincidence probability is found by first 
projecting the two-photon amplitude onto the one fiber, 
and then calculating the overlap with the other fiber in 
the same way as for single coupling. In this example we 
will search for fj.^s and make a conditional measurement 
on the signal, defined by the following operator 

M, = \G^^){G^^\. (21) 

Due to the measurement, the derivation of pi will be 
slightly different here, and we need to take a few steps 
back and reformulate the two-photon density matrix p^j 
as a coherent sum of amplitudes with respect to Acp, in- 
stead of as a incoherent trace operation in Eq. © . The 
density matrix is now written 

pi: = Y.Y.\^tr"'")(^^^'"\- (22) 

Using the measurement operator Mg, the two-photon 
density matrix after the projection becomes 

^'^'1'' TV(Ms ® liP^jMs ® li) ■ ^ ^ 

The reduced density matrix is readily found by tracing 
over the partner, p'r^^ = Trs(pg;|g), which leaves only a 
trace over frequency, pi\s — J2nPi"s- conditional co- 
incidence is now defined in the same way as for single 
coupling; we can replace 7 by ^i|s in Eq. H19|) . still using 
Eq. ifTUI) to find the eigenvalues A„ and eigenmodcs |Cn) 
of pi\s- We have, 

Mi|s = E^»l(^"|Goi])l', (24) 
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where \Gqq) is the fiber-matched mode of the idler. The 
parameter ^s|j follows accordingly, as well as the formal 
optimization: 

/^°P* =max /^(^p,4,i), (25a) 

Ss ,i 

C°P* = argmax /i(ep,Cs,i)- (25b) 
Pair coupling 

Finally, the pair- coupling efficiency 7c is defined as 
the probability to find both photons of a pair in the re- 
spective fiber. This measm^e tells what fraction of the 
pairs enters the fibers compared to the total amount of 
pairs that are generated within the frequency bandwidth 
window. The pair-coupling can be derived from the sin- 
gle coupling and conditional coincidence using effectively 
Baycs's rule, see Fig. |3| 

7c = Mi|s7s = /^s|i7i- (26) 

The alternative is to calculate the coupling via 7c = 
Tr(Ms (g) Mi psi), but this requires the calculation of psi, 
which is computationally more demanding. When com- 
puting ^i|s and 7s via Eq. H2()|) . using Eq. H24|) and 
Eq. ((T^ . the ket is sufficient, because we can simphfy 
the trace-operation of Eq. ||7J), and also the projection of 
Eq. (|23|l . to work in ket-space before the trace over fre- 
quency; pf - Tr^d = ^,n,,5;,, |0^^)(0["^ We 

could also think of rewriting Tr(Ms (8) Mj psi) using two- 
photon kets in the same way, but as psi generally becomes 
a mixture after tracing over frequency this is not an op- 
tion. To compute 7c before the frequency trace is also 
not an option numerically, as the trace over frequency 
involves a for-loop and optimization performed within it 
will reduce efficiency heavily. 

The measure 7c should be compared to 
V = 7c/V7s7i = VMs|iMi|s, which is basically 7c nor- 
malized to 7s and 7i, that have been used by some 
authors [23, |21| . The parameter 77 is useful as a type of 
measure of correlation that tells how well the focusing 
system has been set up to couple the modes of the idler 
emission to the same as those conditioned by the signal 
emission, or vice versa, depending on which of the two 
possess the smaller single-coupling efficiency. We intend 
to simply plot 7c as this compares directly to 7s and 
7i in terms of achievable photon rates; in principle, 7c 
could be low while 77 is high. 

III. NUMERICAL PREDICTIONS 

All results in this section are for the case of a PPKTP 
crystal with the poling period A = 2tt/K = 9.6 pm oper- 
ating at perfect quasi-phase matching; the pump at 532 
nm creates emission at 810 nm and 1550 nm in the abso- 
lute forward direction. The temperature T = 111 °C, 



which affects the fc- vector magnitudes, is chosen such 
that kp = ks + ki + K, see Ref. p^ . 

The numerical calculations are implemented in Mat- 
lab using Eq. Q)-l|21l- AH refractive indices arc deter- 
mined by the ScUmcicr equations |23, [s^j setting the 
wavelength and temperature dependence of the k-vector 
magnitudes. The resolution Ng of the discrete angular 
spectral amplitude representation in the polar degree are 
a few hundred points and varies between 1-100 pradians, 
with the higher resolution for short crystals and strong 
focusing (wide-spread emission) and the lower resolution 
for long crystals and weak focusing (narrow emission). 
The needed azimuthal angle resolution Nip is found to be 
> Ng/5, and the frequency resolution varies between 
a few points for short crystals to a few hundred points for 
long crystals where the spectrally induced contribution to 
spatial multimode is larger. To spare the computer from 
unnecessary workload we observe that the two-photon 
density-matrix in Eq. ((BJ (scaling as Ng number of points 
in size) is always pure and can be fully represented by its 
amplitude vector alone (scaling as Ng), for all of the cal- 
culations. 



A. Single coupling 

As said earlier, according to our definition the single- 
coupling efficiency depends on the emission bandwidth 
filter that is being used. This is because of the fact that 
many of the different frequencies created in the SPDC 
process will not couple into a single-mode fiber. Looking 
at a single frequency of the emission, the angular spec- 
trum of the emission will be described by a single sinc- 
function for each of the plane waves of the pump, see 
Eq. As will be argued in the next subsection, most 
of these sine- functions will overlap nearly perfectly at op- 
timal pump-focusing such that the emission is strongly 
spatially coherent and define almost a single-mode that 
will couple well into a single-mode fiber. If the pump- 
focusing is too weak it will create transverse multimode 
emission, as the many sinc-functions are then distributed 
along the transverse position of the pump beam and do 
not coincide. If the pump is instead focused too strongly 
the effect is the same, except that the multimode now 
originates from longitudinal position, also providing bad 
coupling. This is the general picture using the window 
of a single emission frequency. 

If we look at a wide spectrum of the emission, each of 
the different frequencies can be seen as composed by a 
set of sinc-functions, each set in a different direction, and 
with every sine in a set coming from one plane wave in 
the decomposition of the pump. For long crystals, when 
the width of the sinc-functions narrows down, the differ- 
ent sets of sinc-functions will no longer overlap. Within 
each set the sinc-functions are spatially well overlapping, 
thus defining a coherent single-mode, but as the sets do 
not overlap the emission will become spectrally multi- 
mode similar to above, also resulting in spatial multi- 
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FIG. 5: (Color online) The single coupling of the idler 7°''', 
plotted for a narrow enough filter bandwidth, AAnarrow, which 
shows that about 95% of the emission can be coupled into a 
single-mode fiber at optimal focusing. The solid line shows 
the pump-focusing parameter , and the dashed-dotted lines 
show the focusing of the idler's fiber-matched mode C"'''- For 
each data sample the idler focusing has been optimized for 
maximum coupling using Eq. 1201 . 



mode. This again provides poor coupling efficiencies. 
However, coupling into fibers automatically does some 
spatial filtering as it selects only the coherent part of 
the emission defining a single-mode, i.e. sinc-functions 
largely overlapping, and thereby it also does some fre- 
quency filtering. Altogether, this motivates why we have 
looked at only a single frequency of the emission for the 
results of the numerical calculations of the single coupling 
efficiencies shown in Fig.|5H7| We will refer to this case by 
saying that we have a " narrow enough" filter bandwidth, 
AAnarrow, which maintains a single-mode at optimal fo- 
cusing of the pump and the signal and idler fibers, i.e. the 
bandwidth is narrow enough that the different sinc-sets, 
corresponding to different frequencies, within the band- 
width overlap (are coherent). Frequency filtering effects, 
as those just described, are left to the next section. 

Figurc|Slshows the single-coupling efficiency of the idler 
7i plotted against the crystal length L and the focusing 
of the pump-beam, via its waist wop- For each sample in 
the plot, the idler fiber focusing has been optimized us- 
ing Eq. H20() to find the maximum coupling 7^ As seen, 
there is always the same maximal coupling to be found 
for any length of the crystal by changing the pump-beam 
waist radius accordingly. The straight lines show that the 
focusing parameters of both the pump and the idler 
fiber focusing are constant, which means that the ge- 
ometry of the beam profile and the crystal edges should 
stay fixed for different lengths of the crystal for optimal 
focusing. The said graph would look nearly the same for 
the signal emission, and, taking a different view of the re- 
sults, Fig.Elclearly shows the importance of choosing the 
right combination of focusing for the pump and for the 




FIG. 6: (Color online) The single coupling of the signal 7s, 
plotted for a narrow enough filter bandwidth, AAnarrow, which 
reaches a maximal 98% at optimal focusing, = 1.7 and 
= 2.3. 
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FIG. 7: The single couplings, 7°'' and 7°'' , reaches a maxi- 
mum at ^p = 0.9 for the idler, and at = 1.7 for the signal, 
which corresponds to = 2.4 and 4°p' = 2.3. The line 
representing the signal in this graph is essentially a plot of 
the ridge of the surface in Fig. |S| 



fibers. Interestingly, we observe that as long as the fiber 
focusing is matched to the pump focusing, for any given 
length of the crystal, then the coupling efficiency will 
reach > 45% irrespectively of the pump focusing. This 
fact may very well explain the relatively high efficiency 
nevertheless achived in many fiber-based SPDC-setups 
for which the experimentalist perhaps have not worried 
about changing the pump's focusing, but rather solely 
the fiber coupling. 

Figure |7| shows both the signal and idler coupling in 
a graph that is parametrized by the pump focusing. In 
each case the optimal fiber focusing is found, and plotted 
along the horizontal axis. In this asymmetrical configu- 
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FIG. 8: (Color online) The conditional coincidence plot- 
ted versus the focusing of the pump and the focusing of the 
signal's fiber-matched mode ^s. For each sample in the graph 
the focusing of the idler (5°''* — solid lines) is optimized to 
find the maximum (up to 100%), using Eq. H25^ with a 
narrow signal filter, AAnarrow, and no idler filter. 
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FIG. 9: (Color online) The pair coupling 7c — ^ii|s7s at a 
pump focusing of = 1.3, which is trade-off between what 
is optimal for the signal (^p = 1.7) and the idler (^p = 0.9) 
individually. At optimal focusing, = 2.0 and = 2.3, 
the maximum 7c is about 97%, using a narrow signal filter, 
AAnarrow, aud uo idler filter. 



ration it leads to a maximal 7°p* = 98% when optimiz- 
ing the focusing for the 810 nm emission (^p ~ 1.7 and 
4°P* = 2.3), and 7°^* = 93% for the 1550 nm emission 
(^p = 0.9 and = 2.4). The optimal focusing of the 
pump depends on the amount of non-degeneracy for each 
of the wavelengths, e.g., for the degenerate case (1064 
nm) the optimal focusing is ^p = 1.4 and = 2.3. It 
should be noted that, in general, the found optimal fo- 
cusing parameters do not correspond to a match of the 
beam- waist sizes jl9j . but rather to an equal geometry. 
However, a matching of the waists are within the same 
order of magnitude comparable to using optimal focusing 
parameters. 



B. Coincidence and pair coupling 

For any focusing of the pump-beam, the fundamen- 
tal modes of the signal and idler emission will be highly 
correlated, meaning that, e.g., a signal photon that en- 
ters its fiber will have its idler partner entering the other 
fiber, provided correct fiber focusing. At optimal focus- 
ing of the pump-beam, this correlation is always high if 
the partner beam is focused optimally, independent of the 
focusing of the beam that we condition upon. In other 
words, at optimal focusing of the pump-beam the con- 
ditional coincidence jUj|s, i.e., the probability of having 
the idler photon in the fiber given that the signal photon 
is in the fiber, will be mainly set only by its single cou- 
pling probability 71, which is always at a high value at 
optimal focusing due to the emission being mostly single- 
mode, see Fig. |H1 In contrast, because of the multimode 
character of the emission at other pump-beam focusing 



settings than optimal, a high conditional coincidence can, 
in that case, only be attained near optimal focusing for 
both the signal and idler fibers. Each sample in the plot 
has been generated using Eq. H25|l with a narrow filter, 
AAnarrow, at the signal side, as defined earlier, and with- 
out a filter at the idler side, when finding the maximum 
fj.°^J' that corresponds to optimal focusing of the idler. 

As can be deduced from the graph, the condi- 
tional coincidence is always very high, reaching 100% for 
most weaker focusing conditions. When instead using an 
idler frequency filter that is matched to the signal filter, 
then /ii|s will be bounded above by 71%, assuming Gaus- 
sian shaped filters on both sides. This limitation follows 
from the fact that while the signal photon of a given 
pair may very well be transmitted through its filter, the 
idler may not. Using Eq. Q), the maximum number can 
be easily derived from the normalized overlap integral 
/|A,(e)|2|^,(e)|2de//|As(e)pde = l/^/2, for which we 
note that the result is independent of the bandwidth. 

Additional qualitative results on the optimal joint fo- 
cusing can be found by turning to the pair coupling effi- 
ciency 7c. As opposed to /Xi|s, this measure relates to the 
total amount of pairs that is generated, and not only to 
those conditioned upon. As shown in Fig. |51 for optimal 
pump-beam focusing, there is a maximal value of about 
97% for 7c at = 2.0 and = 2.3. Note that, since the 
optimal pump-beam focusing varies for each of the beams 
for a non-degenerate wavelength case (^p = 1.7 for signal 
and ^p = 0.9 for idler), we had to find a compromise us- 
ing ^p = 1.3. This graph is again plotted using a narrow 
filter at the signal and no filter at the idler. Equation 1261) 
tells us that for matched filters, 7c will also be limited to 
71%, as long as 7s = 1 which is achievable with narrow 
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FIG. 10: The fiber coupled bandwidth is oc 1/L for a wide 
enough spectral filter AAwide, see text, which can be said to 
be the case for the solid line of AA — 25 nm for all crystal 
lengths defined by the plot. In the limit of no filter at all, 
the graph corresponds to the single-mode bandwidth A Asm, 
see Eq. 1281 . The graph shows the result for the signal emis- 
sion (810 nm) at optimal focusing conditions, — 1.7 and 
= 2.4, and the legend shows what filter bandwidth AA was 
used for each line. 
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FIG. 11: The fiber photon flux is oc y/L for a wide enough 
filter AAwido, and oc for a narrow enough filter AAnarrow. 

The filter is defined as narrow or wide in relation to the 
natural single-mode bandwidth A Asm. For the solid line of 
AA = 25 nm the case has been reached where AA = AAwidc > 
AAsM. The graph shows the result for the signal emission (810 
nm) at optimal focusing conditions, = 1.7 and ~ 2.4, and 
the legend shows what filter bandwidth A A was used for each 
line. 



filters. In general, both the conditional coincidence and 
the pair coupling decrease for wide bandwidths; fii\s in 
such case being bounded above by 100% and 7c bounded 
above by the value of 7s. 

In terms of sources of heralded single photons, these 
results imply that almost perfect correlation can be 
achieved by careful focusing and by having no limiting 
interference filter on the triggered photon side; leaving 
such sources limited entirely by the transmission imper- 
fections of lenses and filters, and by detector efficiencies. 



C. Photon-rate and bandwidth 

In this subsection we will look at the achievable photon 
fluxes in free-space and in single-mode fibers and its de- 
pendence on the crystal length. As we will argue, and we 
have shown numerically, this dependence will in turn de- 
pend on the chosen frequency filter. Our arguments will 
follow a series of steps, where the later steps include the 
effects of spatial and spectral filtering. The final results 
are found in Fig. ^| and Fig. ITTI 

As a first step, imagine the pump beam to be a single 
plane wave that is perfectly phase-matched for a single 
frequency of the signal and the idler along the z-axis, 
called here the forward direction. In this case, by look- 
ing at the two-photon amplitude Eq. (^), we see that the 
height of the sinc-function, which describes the angular 
spectrum, is oc L, corresponding to an dependence 
for the intensity. (One should imagine two-dimensional, 
"mexican-hat-like" , sinc-functions.) The width of the 
sine will shrink cx 1/L, such that the fiux will increase 



oc L. This argument is still valid considering the spatial 
transverse multimode emission created by such a plane 
wave pump, discussed earlier. 

As a second step, consider a focused pump being com- 
posed of many differently directed plane waves. In this 
case, still looking at the same single frequency emit- 
ted, each such plane wave will phase-match a little less 
strongly than the one in the absolute forward direction. 
We will have a collection of sinc-functions being added 
together, each originating from a different plane pump 
wave, and numerical calculations show that the combined 
total width, or envelope, of these sinc-functions will de- 
crease for longer crystals, thus adding to the previous 
result a factor 1 / -\/Z, with the flux now becoming oc VX- 

The third step includes the observation that the energy 
of the pump beam is concentrated to the plane wave in 
the forward direction for longer crystals at optimal focus- 
ing. Equation ^ shows that the intensity will be oc uipp, 
because, at optimal focusing we have zn = L/^p, where 
Zfi is given by Eq. H16() . and thus Wq^ oc L. The total 
fiux is now oc 

As a last step we include filtering. In the previous 
steps we looked at a single frequency of the emission, 
which means that the bandwidth was narrow enough for 
the emission to be a single- mode (at optimal focusing). 
For narrow enough bandwidths we therefore get a flux 

P oc L\/ZAA„arrow, (27) 

which is valid both in free-space and in flbcr. As an 
effect of the phase-matching conditions there will be a 
tight connection between the spectral and spatial modes. 
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as we described in Section IIII Al for frequency filtering. 
In terms of fiber- coupling this means that when the fiber 
spatially filters the emission it will also effectively do fre- 
quency filtering. The bandwidth of the signal emission 
(810 nm) coupled into single-mode fibers (using no sep- 
arate frequency filter) is given by 

AAsM = B/L, (28) 

where the value B = 1.23 x 10^^^ [m^] is found for 
PPKTP when both the pump and fiber are focused opti- 
mally, sec Fig.^] We will refer to this bandwidth as the 
single-mode bandwidth. It will also determine how nar- 
row the bandwidth of a filter (AAnanow < AAsm) need 
to be for any given length of the crystal to be considered 
narrow. The photon flux in the fiber will be 

P cx lVIAAsm = VI, (29) 

for any filter A A > A Asm- In Fig. ^2 we have plotted the 
flux for different filters, AAnanow < A Asm < AAwidc- For 
filter bandwidths that arc "wide enough," AAwidc, the 
free-space emission will be multimode even at optimal 
pump focusing, and the free-space photon flux becomes 

Pcx \/Z5(AAwidc), (30) 

where g is some unknown and non-trivial function deter- 
mined by the properties of the crystal material via the 
Sellmeier equations. 

These results clearly show that it is advantageous to 
have long crystals as the photon-rate will always mono- 
tonically increase even when coupling the emission into 
single-mode fibers. As an effect, we can keep the pump 
power low. promoting the use of a compact and cheap 
laser. This requires that we change the focusing of both 
the pump and the fibers ^s,i to the optimal for some 
length L. Additionally, longer crystals give narrower 
bandwidth, which is very advantageous in many appli- 
cations of entangled photons. For example, in time- 
multiplexed schemes it is crucial that the photon pack- 
ets keep their widths in the fibers and do not broaden 
due to chromatic dispersion, and the broadening can be 
limited by having a narrow bandwidth. Another way 
of reducing the effect of broadening is by introducing 
negative dispersion using an appropriatly designed fiber 
Bragg grating. In general these have to be custom man- 
ufactured for broad bandwidths, but for telecom band- 
widths, 30-80 GHz, (in the C-band, between 1525-1562 
nm) these are standard off-the-shelf items, and corre- 
sponds to wavelength bandwidths of about 0.25-0.65 nm 
at 1550 nm. Wc can see from Eq. (^5)1 that 70 — 180 mm 
long crystals are needed, taking into account the conver- 
sion factor between signal and idler bandwidths [AAj = 
(Aoi/Aop — l)^AAs ~ 3.66 x AAg]. Narrow bandwidth 
can of course be obtained by the use of spectral filters, 
however, our results show that it is better in terms of 
photon-rates to use long crystals to achieve small band- 
widths rather than to strongly filter the emission of a 



Beam quality factor 
30 l 

20- = :;:;:;:;:;;;::: = :::;:;:;;;;:;::::;::;::;:;=^^ 




^ 0.01 0.1 1 10 100 1000 

% 

FIG. 12: The beam quality factor of the idler plotted 
against the pump beam focusing ^p. The smallest value, 
= 1.4, is found for = 0.9 



short crystal. (This is in contrast to what is claimed by 
Lee et al. in Ref. |33j. for birefringent phase-matching 
and intersecting cones.) Furthermore, with narrow band- 
width follows also long coherence length of the photons 
which is highly desirable when working with intcrfcrom- 
ctry as is commonly done when using time-multiplexing 
analyzers to code and decode qubits. 



D. M"^ and coupling 

In this subsection wc will present the numerical predic- 
tions of the emission mode in terms of the beam quality 
factor for different focusing conditions. We will also 
elaborate on the connection between the beam quality 
factor and the coupling efficiency 

Figure ^] shows the beam quality factor M? plotted 
against the focusing of the pump for a narrow enough 
frequency bandwidth of the idler emission (AAnarrow ^ 
AAsm)- There is a clear optimal focusing, where the 
emission reaches close to single-mode, M? ~ 1.4, at a 
focusing of = 0.9. These results are valid for any 
length of the crystal, compare to Fig. [S] A low value 
of means that the light is close to a single-mode, 
and thus possible to couple well into a single-mode fiber. 
For bandwidths larger than the single-mode bandwidth 
AAwidc 2> AAsm, the light will become spatially mul- 
timode and the coupling efficiency will decrease accord- 
ingly 

Figure [TSI shows the relation between the coupling ef- 
ficiency 7i and the M^, as the focusing of the pump 
is varied. The correspondence is clear, and we can see 
that different values can provide the same coupling 
efficiency. This is so because the coupling efficiency is 
only determined by how much of the emission is in the 
fundamental mode. What determines the Af^ is the dis- 
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FIG. 14: The experimental setup used to create polarization 
entangled photon pairs, and to verify numerical results. PBS: 
polarizing beam splitter; HWP: half wave plate; QWP: quar- 
ter wave plate; SWP: short-pass filter; BP: band-pass filter; 
SMF: single-mode fiber; 77: detection efficiency. 



FIG. 13: The single coupling 71 versus the of the idler, 
using the same data as in Fig. 1121 and Fig. [7| The graph is 
parametrized by the pump beam focusing and illustrate how 
a low is connected with a large 71. 



tribution of the light between the higher order modes, 
and this can differ from one case to another, even with 
the same amount contributing to the fundamental mode. 
In general, as we have said, too weak focusing will pro- 
vide spatial transverse multimode, and too strong focus- 
ing will provide spatial longitudinal multimode. It can be 
deduced from Fig. I13l that longitudinal multimode, orig- 
inating from too strong focusing, creates emission with 
relatively higher contribution to the fundamental mode 
for the same A''P value. 



IV. EXPERIMENTAL RESULTS 

To verify some of the numerical results we compared 
with experiments. We have measured the beam quality 
factor, the bandwidth in the fiber, and the coupling ef- 
ficiencies for different focusing conditions of the pump. 
The experimental setup is shown in Fig. 1141 As a pump 
we use a frequency doubled YAG laser emitting approx- 
imately 60 mW in the TEMqo mode at 532 nm. Its 
Mp-value was measured to 1.06. After a band-pass fil- 
ter (BP532), which removes any remaining infrared light, 
we "clean up" the polarization using a polarizing beam 
splitter (PBS). The polarization is controlled by a half 
wave plate (HWP) and a quarter wave plate (QWP) in 
front of the crystal. The pump-beam is focused onto the 
crystal using a achromatic doublet lens (/p = 50 mm) 
which introduces a minimal amount of aberrations not 
to destroy the low value. The QWP is set to undo 
any polarization elliptisation effects caused by the lens, 
and fluorescence caused by the same lens is removed by 
a Schott filter (KG5). 

The next component is the crystal. This is a peri- 
odically poled, bulk 4.5 mm long KTP crystal, with a 
poling period of A = 9.6 fim, which will colinearly cre- 



ate a signal at 810 nm and an idler at 1550 nm when 
heated in an oven to a temperature T sa 100°. When the 
setup is used to create polarization entanglement, two 
crystals are present, one oriented for V and one for H, 
and the polarization of the pump is set to 45°. By cou- 
pling the emission from both crystals into single-mode 
fibers we cannot even in principle determine which crys- 
tal the photons came from, except by their polarization 
degree of freedom, and therefore the signal and idler will 
interfere in the diagonal basis and get entangled in polar- 
ization. This principle was first demonstrated by Kwiat 
et al. in Ref . [23 ■ Our first results was presented in Ref. 
[2^ . and the latest results, overcoming some problems 
of crystal dispersion and using optimal focusing, will be 
found in Ref. 0. 

After the crystal, we block the pump light by a 532 
nm band-stop filter, and the signal and idler emission is 
focused by achromatic doublet lenses. The rather small 
F-number [F = //!?, where / is the focal length and D 
is the beam diameter) of the emitted light {F < 40 for 
/p ~ 50 mm and F < 9 for /p = 12 mm) requires good 
quality lenses not to increase the factor. The lenses 
we use are all aberration free down to F « 6 — 11, and 
are also quite insensitive to an offset in the alignment of 
the optical axis. 

To determine the coupling efficiencies and bandwidths, 
the complete setup of Fig. [HI was used. To separate the 
810 nm and 1550 nm emission we used a dichroic mir- 
ror made for a 45° angle of incidence. The first lens 
(/si = 30 mm) is common to both signal and idler and 
its task is to refocus the beams somewhere near the 
dichroic mirror. The next two lenses (/g = 60 mm and 
/i = 40 mm) coUimate each beam, and they are focused 
into the fiber-tips (with the mode field diameters being 
MFDsio = 5.5 fim and MFD1550 = 10.4 ^m) using as- 
pherical lenses with / = 11 mm. In front of the fiber 
couplers we have first Schott filters (RG715) to block 
any remaining pump light, and then interference filters 
of 2 nm and 10 nm at the 810 nm and 1550 nm side 
respectively (BP). The detectors used were a Si-based 
APD (PcrkinElmer SPCM-AQR-14) for 810 nm and a 
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FIG. 15: The experimentally observed beam quality factor, 
M^, for the idler beam at different sizes of the pump beam 
waist radius luop- The lowest value of the is 2.8 at a 14 ^m 
pump waist. 



the pump lens: 12 mm. 30 mm, 50 mm, 75 mm, 100 mm, 
and 150 mm, each being placed at a distance that set the 
focus in the center of the crystal. The result, which is 
shown in Fig. 1151 agrees fairly well with the numerical 
predictions. The shortest focal length lens, 12 mm, gave 
a somewhat higher M^, which can be explained by the 
fact that this was the only singlet lens used, probably 
adding some aberrations, while the others where achro- 
matic doublets. The lowest value, M?^ — 2.8, was found 
with the 50 mm lens giving a 14 /im pump waist radius 
wpp inside the crystal, corresponding to = 2 for the 
4.5 mm long V-crystal (for later reference we observe that 
= 1.3 for L = 3 mm agrees a little bit better with nu- 
merical results). Note that the values are slightly 
higher here compared to Fig. ^1 This can be explained 
by the non-perfect phase-matching in the experimental 
case, resulting from cither too low crystal temperature, 
uncertainty in the true value of the poling period (possi- 
bly deviating somewhat from its specification), or both. 



homemade InGaAs-APD (Epitaxx) module for 1550 nm. 

When determining the beam quality factor, M^, we 
used only a single crystal oriented to create vertical (V) 
polarized light, and the complete setup of Fig.lTHwas also 
not used. Instead, we focused the idler emission directly 
using a lens of focal length j\ = 75 mm placed at a dis- 
tance of 75 mm from the V-crystal to collimate the beam. 
At the additional distance of 470 mm we placed another 
lens with focal length /i = 150 mm that refocused the 
beam again, so that we could take measurements of the 
beam profile around its waist. 



A. measurements, results 

To obtain the results of Fig. ^| we first took images 
of the refocused idler beam in the x-y plane using an 
InGaAs-detcctor camera from Indigo Systems, model Al- 
pha NIR. Several images were acquired for different po- 
sitions along the z-axis around the waist, and we then 
integrated the resulting 2-dimensional surface over one 
axis to create an intensity profile for the remaining axis. 
Because of the detector noise we could not use the stan- 
dard deviation method to find the beam radius, defined 
by the 1 /e^ level. Instead we matched a Gaussian shaped 
function to the intensity profile to find the width. This is 
accurate enough for mode-shapes that arc close to Gaus- 
sian, which is the case for low AP values. To limit the 
impact of the noise we applied a function that assigned 
greater weight to the center- values of the intensity pro- 
file. The widths of the beam for each z-axis position 
were then set together to find the beam profile of the 
emission, and its factor was determined by fitting to 
the standard Gaussian-beam function, Eq. I|17|l . We now 
repeated the procedure for different focal lengths, /p, of 



B. Coupling efficiencies, results 

The experimental data for the coupling efficiencies 
were obtained with the source producing polarization- 
entanglement using two crystals. For this reason we ex- 
pect the values to be a bit lower than predicted as we 
needed to focus the fiber-matched modes for both the H 
and the V crystal at the same time. We also have this 
problem with the pump beam, and we aimed at placing 
the focus at the intersecting faces of the two crystals for 
both the pump and the fiber. As already mentioned, the 
temperature of the crystal used in the experiment was set 
lower than required for absolute perfect phase-matching 
at 810 nm and 1550 nm. This was because we observed 
higher photon fluxes at this setting. Contradictory as 
it may seem, the explanation is that the peak of the 
emission spectrum is not symmetrically centered around 
the above wavelengths, but rather towards 810 — a and 
1550+/3, including a long tail representing the emission at 
larger angles. As our filters are centered for 810 nm and 
1550 nm, the peaks of the emission can be moved to line 
up with these by changing the temperature, and thus the 
phase-matching, which will give somewhat higher fluxes 
although the coupling efficiencies will decrease accord- 
ing to our definitions. In addition to having a slightly 
wrong poling period these effects degrades the efficien- 
cies, which we could verify numerically and which is sup- 
ported by comparing Fig. 1151 and Fig. ^1 The obtained 
results for the single coupling efficiencies were 7s = 32% 
and 7i — 79%, for the conditional coincidence /ii|s — 34%, 
and for the pair coupling 7c = 11%, when focusing ac- 
cording to = 2.1, £_s = 3.2, and = 2.5 (as decided by 
available lenses, and assuming L = 4.5 mm). For these 
numbers we have compensated for the 35% transmission 
of the 1550 nm filter, and the 85% transmission of the 
810 nm filter. The singles photon rate in the signal fiber 
was 2.3 Mcps (10^ counts/sec) and in the idler fiber 2.4 
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Mcps. The total generated rate of photons before fiber 
couphng was estimated at 8.6 Mcps and the coincidence 
rate in the fibers was 274 kcps, (see Ref. |23|)- 



C. Bandwidth, results 

We have used a spectrograph (SpectraPro 500i, ARC) 
to measure the bandwidth of the signal emission using 
the single-mode fiber without a filter. The bandwidth 
was 4 nm for the V-crystal and 6 nm for the H-crystal. 
Figure El suggests that the effective length of the crystal 
being poled must be 3 mm and 2 mm respectively. Also, 
from Fig. ^] for the 2 nm filter, we can deduce that 
the 2 mm crystal should give roughly 55% of the photon 
rate of that of the 3 mm one. Experimental agreement is 
good, as we saw the H-crystal giving half the rate of the 
V-crystal (with no compensation done by balancing the 
fiber coupling or rotating the pump polarization). Re- 
ferring again to Fig. El using the effective crystal length, 
the best pump beam focusing parameter is modified to 

= 1.3 for L = 3 mm (V-crystal) which agrees roughly 
with the value of optimal focusing, = 0.9. 



V. CONCLUDING DISCUSSION 



liminary numerical calculations then show increased 
values and decreased coupling efficiencies. However, due 
to the characteristics of the convolution, it seems we can 
retain the good results of CW even for pulsed operation 
by using very long crystals, as this will bring back the 
delta-function at the limit of infinitely long crystals. For 
this discussion wc have not yet worried about any disper- 
sion effects that might come with long crystals and short 
pump pulses. 
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APPENDIX A: THE TWO-PHOTON 
FREQUENCY AND ANGULAR SPECTRAL 
AMPLITUDE 



In summary, precise focusing of the pump-beam 
and the fiber-matched modes can significantly increase 
the coupling and coincidence efficiencies of quasi-phase 
matched SPDC-sources, which is important for applica- 
tions needing highly correlated pairs of single photons to 
propagate in fibers. We have shown how the beam qual- 
ity factor of the emission changes with the focusing of 
the pump. At optimal focusing the emission is mostly 
created in a spatial single-mode, which couples well into 
single-mode fibers, and by maintaining a fixed geome- 
try of the beam profile in relation to different lengths of 
the crystal this stays true for all lengths. We have also 
shown how the photon flux depends on the crystal length 
for different frequency filters, the conclusion being that 
longer crystals produce more photons per unit time at a 
smaller bandwidth. 

In all of the calculations we have assumed a monochro- 
matic (CW) pump laser. Looking for a possible extension 
to pulsed operation we observe that the interaction time, 
r, in Eq. ljA18|) for a CW laser is set by the coherence 
time of the pump alone, and as T is infinite it trans- 
forms into a delta- function of frequency in Eq. I|A20(I . 
Using pulsed light, the integral exp(— iAwt) should 
be replaced by /i(t) exp(— i Awt), where h{t) is the 
convolution, h{t) = hc{t) * hL{t), between the form of 
the temporal wave-packet of the pump, hc(t), and the 
form of the crystal along the 0-axis, h^lt). We observe 
that when hc(t) is narrow, like for pulsed operation, the 
transform of h(t) will instead become a sine- function, 
specifying an inexact energy-matching condition. Pre- 



The evolution of the number state vector is given by 

to+T 



(Al) 



where |'0oo) is the state at time to, T is the time of inter- 
action, and H{t) is the Hamiltonian 




H{t) = j x'^^^E^+^Et^E[-^d\ + H.: 



(A2) 



There are three interacting fields in the crystal's volume 
V ignoring all higher-order terms (n > 3) of the nonlin- 
earity x^^^ ■ AH three fields have the same polarization 
{ZZZ): 

eI+^ = J2 Ap(sp)e*('=p^p-''-'^p*+"^p\ (A3a) 

Et^ = Jd(bsJdLo,A{Lu,) ^ e-*('==«=-'-'^=*+*=)at(^„ s,), 

(A3b) 



E: 



(-)_ 



(A3c) 
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The field of the pump is classical and monochromatic so 
that we can replace £^p~'''' by Ep'^\ The plus-sign de- 
notes conjugation, i.e. annihilation (+) or creation (-) 
of the state. In all the calculations we use the notation 
k = ks, where s is the unit length vector of k. The an- 
gular amplitude spectrum Ap{sp) takes into account the 
focusing of the pump. For signal and idler, we sum over 
both frequency and angular modes, where a(w, s) is the 
field operator, and A{iu) is the frequency amplitude of 
a Gaussian shaped detector filter having the bandwidth 
AA (FWHM) and center wavelength (all wavelengths 
in vacuum). Via the relation uj = 2TTcn^/X its form is 
given by 

A(a;;A)==e-2iog(2)(A-AjVAA^, (^4) 

Each signal and idler photon is created with a random 
phase, 0s and 4>i respectively, which we also need to sum 
over. The only nonzero solution is completely correlated 
phases as will be shown later. The phase of the pump (j)p 
is constant but arbitrary. 

For periodically poled materials, the nonlinearity x'^' 
has sharp boundaries, and later on in the calculations 
it will facilitate to make an expansion of x^^^ into its 
Fourier-series components 



X 



(2) 



X2 fir) = X2 XI 



-imK-r 



(A5) 



Tn=0 



and then do a sinusoidal approximation using the first 
term, 

X^'^ = X2 fie-''' '-, (A6) 

where K = 2tt/A e^, and A is the grating period. Ap- 
pendix^treats the case of a M+1 term series expansion. 
From Eq. (jAl|) the number state becomes 

IV") = iV'oo) + JJdujsdcoi'^'^ S{ujs, uji, Ss, Si)dlal\il}oo) 



(A7) 



where G2 is the unnormalized amplitude for the two- 
photon number state. 



G2 = (V'lilV') = / duJsduji S{uJs, LOj, Ss, Sj), 



such that for to — 0, 

T 

1 



I dt H{t) = G2 alal - H.c. 



(A8) 



(A9) 



Our goal now is to arrive at an expression for the am- 
plitude S which will also enter in the state of frequency 
and angular spectrum of the form 



(AlO) 



We start by inserting Eq. ljA6p into Eq. (|A2|) and then 
Eq. HA2|I into Eq. which gives 

T 

G2^^Jdt jd\ X2 fie-'''-''El+^Et^El'\ (All) 
V 

By making a substitution of the fields in Eq. (jA3|) into 
Eq. (|A11|1 . and via identification using Eq. I|A8|) we find 
that 

X2 flA{uJ,)A{u!i)J2M^p) 

L/2 oQ 00 
X Jdzjdyjdx e-'^k-(xa^+y<,,+za,^) 

— L/2 -00 -co 

2iT T 

^ ^JJ d0sd0i J dt e-»[('^=+'^.-'^p)* + •P^+'i'^-^'p]^ 


(A12) 

where the volume integral has been expressed in a Carte- 
sian coordinate system (r = xe^ + yey + ze^, see Fig.^, 



L/2 



00 00 



d r = J dz J dy J dx. 

V -L/2 -00 -00 



(A13) 



We have also introduced the phase mismatching vector 



Afe = kgSs + kiSi — fcpSp + K 
+ AkyBy + AkzBz 



(Al4a) 
(A14b) 



In a Cartesian coordinate system the normalized vectors 
s arc represented by 

Ss = Ps e^. + Qs By + rris e^. 
Si = Pi e^. + qi By + mi e^, 



Pp ex + Qp By + nip Bz 



K = K Bz 



(A15) 



where p, g, and m are the normalized components of s 
in each of the three dimensions 113 • 

Because of the rotational symmetry of the emitted 
modes, it is suitable to use a spherical coordinate sys- 
tem {0,(p), for which p = sin9 coaip, q = sin sin t^, and 
m = cos 6. The phase-mismatch vector components then 
become 

Akx = ks sin 9s cos ips + h sin 9i cos ipi — hp sin Op cos Lpp, 
Aky = ks sin 9s sin (ps + ki sin 9i sin ipi — kp sin Op sin ipp, 
Akz = ks cos Os + h cos Oi — kp cos Op + K. (A16) 
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Note that the magnitude of the signal and idler fc- vectors 
implicitly depends on the polar angle 9 according to 



KiOs) = 1/1 



k-m = 1/1 



cos 9s ^ 


2 

+ 


/ sin6's 




I kl 


COS 9\ \ 




( sin^i 




^ kY 



(Al7a) 
(AlTb) 



where fcf , fc^, , and \i( are the constant magnitude 
of the fc- vectors along the crystals Z and Y axis, respec- 
tively (fcp need to be constant and equal to k^ as we 
will soon show). Generally, there is negligible difference 
in refractive indices between the crystal's X and Y axes 
which cancels the dependence on the azimuthal angle ^) 
in the equations above. We therefore use the Y axis as 
the major axis being orthogonal to Z . 

Using spherical coordinates exclusively leads to 

5(tJs,t^i,6's,6'i,'y3s,</2i) = 

7r/2 2ix 

Xi /iA(cjs)^(wi) y sin6'p d9pjdipp Ap{9p,tpp) 


L/2 00 00 



-L/2 -00 —00 

2tt T 

1 

ih 



(A18) 

The angular spectral amplitude Ap of the pump beam in 
Eq. (|A18p is Gaussian shaped for a laser emitting in a 
TEMqo single-mode, and in spherical coordinates it be- 
comes |30l 



p/4 



2tt 



(A19) 



where the beam waist radius WQp of the focused pump- 
beam has entered the calculations. The function is nor- 
malized to represent the same constant power available 
in the beam at different focusing conditions. 

Now we will solve the integrals over space, time, 
and phase in Eq. IjAlSp . In doing so we note that 
there are three spatial integrals of which two are the 
Fourier transforms of unity (dx and dy) and one is the 
transform of a box- function (dz). The transforms turn 
into two (5-functions and a sinc-function respectively. 
The time-integral also turns into a (5-function of the 
three frequencies cOg, loi, and cUp. This is because 
we have a monochromatic pump-beam with infinite 
coherence length, which effectively leads to an infinite 
interaction-time, T — > 00, even for short crystals. The 
two integrals over the random phases (j)^ and (f>i will 
make the amplitude S vanish completely if the phases 



are not fully correlated with each other. Therefore, the 
only nonzero solution is when the two phases add up to 
a constant. S can be complex-valued, thus yielding the 
relation (p^ + 4>i ^ 't'p + C . If we let C = for simplicity, 
we are led to 

S{LOs,i^i,9s,9i,ips,<Pi) = 

7r/2 2tt 

X2 fiA{ujs)A{uJi) j sin6'p d9pjdipp Ap{9p,(pp) 




X S{Akx) 6{Aky) L sine 
47r2 

X -—5{ujs+i^i-i^p)- 
in 



2 



(A20) 



We now have two integrals over 9p and ipp with S- 
f unctions over Ak^ and Aky which in turn depends on 
9p and (pp according to Eq. ljA16|l . The integrals can be 
canceled in a few steps by setting the equalities Ak^ = 
and Aky = 0, and to that end we need to assume that 

fcp which 



fcp is constant for small angles 9p, i.e. fc 



we believe is a fair approximation for pump-light that 
is not extremely focused. By extreme we mean beyond 
the validity of the paraxial approximation. The latter 
equality applied to Eq. ljA16|) gives 



kg sin 9s sin tps + fci sin 0; sin ipi 
fcp sin 9'p 



(A21) 



Equation ljA2ip together with the relation arcsin(a;) = 
arccos(-\/l — a;^) now gives the following expression for 
Afc^ = of Eq. IjAiep (with ipp primed), 

fcg sin 9s cos ips + fci sin 0; cos ipi+ 



-k^sm9'pJl- 



ks sin 9s sin tps + fci sin 9i sin cpi 



fcp sin 9'p 



= 0. 
(A22) 



If we now take the square of Eq. HA22() and solve for 9p 
we get 

9' ~ arcsin \/ P'^ + Q'^ = arccos \/\ — (P^ 



where 



P 



Q 



ks sin sin ips + fci sin 9i sin Lp{ 
/tp 

fcs sin cos (ps + fci sin 9{ cos ip{ 



uz 
ftp 



Q2), 

(A23) 

(A24a) 
(A24b) 



Furthermore, 



p2 + 1 



fcg sin^ 9s + kf sin^ 9i + 2fcsfci sin 9s sin 9i cos(A(/3) 



(A25) 
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where we are allowed to introduce A(/3 = ips — ipi- This 
is a result of the assumption of rotational symmetry and 
will lead to the final state being invariant to a common 
variation in the azimuthal angles for signal, ips, and idler, 
Lpi. As shown here, only the angle-difference is of im- 
portance. Using Eq. l|A23p in the expression for Afc^ of 
Eq. IjAlGp wc have 



APPENDIX B: SERIES EXPANSION OF X^^^ 

The poling structure of periodically poled crystal has 
the approximate form of a square-function along the z- 
axis. In such a case, the M + 1 term series expansion of 
x'^^ become 



Afc; fcs cos 6*3 + fci cos Oi ^ yfl - {P'^ + Q^) + K. 

(A26) 

At this stage the two integrals in Eq. (|A20|) have been 
canceled and the amplitude can be simplified as 



X2 /iA(u;s)A(a;i)Ap(0' 



X L sine 



-Afc^ 

2 ^ 



47r2 

X — — (5(tJs + cji - cjp). 



{A27) 



X^^' = X2 fir) = ^ E i-il::^e-(2™+i)^-, (Bl) 
TT — ' Zm + 1 

m=0 



where K = 27r/A e^. and A is the grating period. In 
the following expression wc have isolated the z-dependent 
part of Eq. (|IT^ : 



L/2 

X2./1 J e" 



(B2) 



-L/2 



One further simplification includes the observation 
that the frequency (5-function can be reduced to unity 
by introducing a common frequency e instead of ajg and 
LUi as defined by uJs = luqs + e, cjj = luqi — e, so that for 
two matched filters the form of the filter amplitude 
becomes squared. Using also Eq. HA23() together with 
Eq. l|IT9)l the expression for the amplitude of the state 
of frequency and angular spectrum finally becomes 



4n X2hL ^^2(^^yp 



X smc 



ih 
L 
2 



(ks cos 0s + h COS 6*1 - k^ - (P2 + Q2) _^ X 

(A28) 



25I and the fcs's and 



where -I- is defined by Eq. 
fci's by Eq. (Tn?!) . 

We now have a final expression for the two-photon am- 
plitude 



G2= /de 



; j I'sm6sd6ssm9id9i JdAif S{e,es,ei, Aip), 

(A29) 



which gives the two-photon state-vector in terms of fre- 
quency and angular spectrum in the form of Eq. ljA10|l 



\A,eAv)=G2\eMm\A^). 



(A30) 



Now, putting the series expansion of x^^-* into the calcu- 
lations of Appendix [K\ the former expression should be 
replaced by 



4X2 

TT 



L/2 



M 



where 



-L/2 



Afcl' 



m=0 



(~l)'"_g-iAfc' 



(™), 



2m +1 



(B3) 



Ak' + 2mK. 



(B4) 



By reversing the order of the sum and the integral in 
Eq. IjBSp we can identify a Fourier transform of box- 
function with an extra phase. The result of the transform 
is a sine, providing thus 



4X2 ^ i~iy 



TT ^ 2m -(- 1 

m=0 



^{Ak',+2mK) 



(B5) 



which is the final expression to replace the sinc-function 
in the state amplitude, Eq. (|A28p . having now M + 1 
terms to approximate the square-shaped poling struc- 
ture. For M = the expression reduces to the sinusoidal 
approximation with fi = A/n. 
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